Reconciling the correlation length for high-spin Heisenberg antiferromagnets 
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We present numerical results for the antiferromagnetic Heisenberg model (AFHM) that definitively 
confirm that chiral perturbation theory, corrected for cutoff effects in the AFHM, leads to a correct 
field-theoretical description of the low-temperature behavior of the spin correlation length for spins 
S > 1/2. With two independent quantum Monte Carlo algorithms and a finite-size-scaling technique, 
we explore correlation lengths up to ^ ~ 10^ lattice spacings a for spins 5 = 1 and 5/2. We show how 
the recent prediction of cutoff effects by P. Hasenfratz is approached for moderate ^/a — 0(100), 
and smoothly connects with other approaches to modeling the AFHM at smaller correlation lengths. 
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In the past decade there has been a resurgence of in- 
terest in the quantum Heisenberg model This interest is 
mainly due to the discovery that the undoped, insulat- 
ing precursors of lamellar high- Tc superconducting cop- 
per oxides are well-described by the spin S — 1/2 two- 
dimensional (2-d) antiferromagnetic Heisenberg model 
(AFHM) on a square lattice. The field theories that de- 
scribe this model at low temperatures share the property 
of asymptotic freedom with the theories that describe ele- 
mentary particles, thus earning a share of attention from 
the high-energy physics community as well. 

The low-temperature physics of the AFHM is dom- 
inated by magnons. The magnon interactions are de- 
scribed by the 2-d classical continuum 0(3) non-linear a- 
model at large correlation lengths. This is an extensively 
studied model in field theory, offering various known ex- 
act results that can be exploited for the prediction of 
the correlation length S,{T) in the Heisenberg model at 
low temperatures. The challenge is to find a proper way 
to connect the parameters of the quantum Heisenberg 
model with the coupling of the cr-model. Several ap- 
proaches to this problem exist. Chakravarty, Halperin, 
and Nelson [I| used renormalization group arguments to 
predict the leading behavior of ^(T). Hasenfratz and Nie- 
dermayer |^,g| utilized analytical results for the c-model 
to refine the prediction, and they used chiral perturba- 
tion theory (CPT) g for the AFHM g to connect the 
parameters of the two models. 

Neutron scattering experiments on S* = 1/2 antiferro- 
magnets such as Sr2Cu02Cl2 generally agree with this 
prediction H. However, higher-spin antiferromagnets 
have remained problematic. Widely different techniques 
- including experiment R~M1 > high-temperature series ex- 
pansion pOl , quantum Monte Carlo simulation for S = 1 
| pd| , and semi-classical approximation |l2| - all showed 
large deviations from the field-theoretical prediction by 
as much as 75% for ^ < 200 lattice spacings, which is the 
regime accessible in experiments on La2Ni04 (5 = 1) 
and Rb2MnF4 (5* = 5/2). These results suggest that a 



serious discrepancy would persist to really large, macro- 
scopic correlation lengths for S > 1/2, which is highly 
unsatisfactory on theoretical grounds. 

In a recent paper ||l^ , Hasenfratz argues that this dis- 
crepancy is due to cutoff effects in the AFHM, which in- 
crease strongly with spin 5*. When large, they can not be 
described anymore with the effective approach of CPT. 
Hasenfratz used spin-wave expansion to calculate the cut- 
off effects, and he showed the proper way to incorporate 
them into the CPT result. 

In this Letter we show via extensive quantum Monte 
Carlo (QMC) calculations that this correction indeed ac- 
counts for the severe spin dependence of the correlation 
length. Our data connect the regime of large and mod- 
erate correlation lengths, where "CPT-|-cutofF' applies, 
with the regime of small correlation lengths where high- 
temperature and semi-classical results apply. The diverse 
approaches are thereby reconciled. We also provide ev- 
idence that by S* = 5/2 the residual deviation at small 
correlation lengths, possibly due to the missing higher- 
order terms in the analytical calculations, has essentially 
reached the classical S' — > cx) limit. 

Consider the 2-d quantum Heisenberg model with 
nearest-neighbor interaction on an L x L lattice with lat- 
tice spacing a and periodic boundary conditions. 



H — J 2_^ Sx ■ Sx+fi , S^ — S{S + 1) 



(1) 



where J > is the antiferromagnetic exchange coupling, 
/t denotes the two primitive translation vectors of the 
unit cell, and Sx is the spin operator at position x. 

The 0(3) symmetry of this Hamilton operator is spon- 
taneously broken at zero temperature, and the model ex- 
hibits long-range antiferromagnetic order in the ground 
state. As a consequence, the model has two massless, 
relativistic Goldstone bosons (called magnons or spin 
waves). However, the Mermin-Wagner-Coleman theo- 
rem [|4| rules out Goldstone bosons in two dimensions 
at nonzero temperature. Instead, the AFHM magnons 



acquire a mass c/£,{T) where c is the spin-wave veloc- 
ity. (We set h, and fc^ to unity.) In fact, the cr-model 
is known to have a non-perturbatively generated mass 
gap, and the leading exponential behavior of the correla- 
tion length in the AFHM is a consequence of asymptotic 
freedom in the cr-model. 

The partition function of the 2-d quantum spin model 
in Eq. (^ can be represented by a path integral jl5| of 
a classical model with an additional ( "time" ) dimension. 
The continuous coordinate X3 of this periodic Euclidean- 
time dimension has extent c/T. When T — > 0, the cor- 
relation length grows exponentially, and becomes much 
larger than the length scale c/T. The system then ap- 
pears dimensionally reduced to a thin slab with two in- 
finite space directions and an extent c/T ^ £,{T) in the 
Euclidean-time direction. This is just a special regular- 
ization of the cr-model. 

Hasenfratz and Niedermayer [g used CPT for the 
AFHM Q, as well as the exact mass gap Q and the 3- 
loop /3-function |]16[ of the cr-model to derive the asymp- 
totic prediction for the spin correlation length in the 
AFHM: 
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The values of spin stiffness ps and spin-wave velocity c 
are not fixed by CPT, but they can be estimated by, for 
example, spin-wave expansion (SWE) ||l^ or fits from 
QMC data [|l|,0. Calculating the 0(T'^) corrections 
in CPT introduces new, unknown parameters. We call 
Eq. (||) the CH2N2 formula after its parents Chakravarty, 
Halperin, Nelson, Hasenfratz, and Niedermayer. 

The QMC data presented in this Letter confirm that 
the discrepancy between CH2N2 and the S > 1/2 AFHM 
correlation length indeed is severe, and it persists to 
macroscopic correlation lengths £,/a ^ 10^. This sit- 
uation is unsatisfactory since a mapping of the AFHM 
onto the cr-model must be valid for smaller correlation 
lengths, too, due to the dimensional reduction described 
above. In particular, this mapping is valid beyond the 
regime of "renormalized classical scaling" pffl near T = 
(where Eq. (0) unquestionably is valid for any S). The 
virulent discrepancy indicates a shortcoming in the tech- 
nique that connects the coupling of the cr-modcl with the 
AFHM parameters ps and T using CPT for large spin. 

In his recent calculation ||l3|, Hasenfratz used bilinear 
spin-wave expansion to modify this connection, taking 
into account cutoff effects in the AFHM. In the present 
study, we also account for a minor refinement of the re- 
sult by Hasenfratz: a part of the quadratic temperature 
dependence, coming from known terms in the spin-wave 
expansion of the coupling of the cr-model, is incorporated. 
The resuhing CH3N2B formula is (with t = T/{2ttPs)) 
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The parameter 7 = 2JS/T brings in the explicit spin 
dependence. In Ref. iQ, exp(— (7(7)) is expressed as 
an integral of familiar spin- wave quantities over the first 
Brillouin zone. The asymptotic T — > (S* fixed) behav- 
ior is C(7 ^ 00) -- 7~^, so the CH2N2 formula, Eq. (||), 
is recovered in this limit. The effect of the aforemen- 
tioned refinement is simply to add the term (27/32)i^ to 
the polynomial. (This term is only a part of the 0{T^) 
correction.) 

In Ref. [gl|, an efficient continuous Euclidean-time 
QMC algorithm ||I|,|l|,||l was used to study the S = \/2 
AFHM correlation length up to 350,000 lattice spac- 
ings. For this purpose, a finite-size-scaling technique, 
developed by Caracciolo et al. ("CEFPS") || for the 
cr-model, was applied to the AFHM finite-volume ^{L) 
data. That study confirmed the validity of the "no cut- 
off effects" CH2N2 formula, Eq. (||), but only for large 
^/a « 10^. By fitting a naive quadratic term at^ in 
the polynomial factor of CH2N2, good agreement was 
achieved down to £,/a « 100 (yielding a = —0.75(5)). 

For S > 1/2 we used two independent QMC al- 
gorithms. One is a higher-spin generalization of the 
S — 1/2 continuous-Euclidean-time loop-cluster algo- 
rithm first described in Ref. |13|. The other one is a 
"traditional" discrete-Euclidean-time loop-cluster algo- 
rithm for arbitrary S, based on a method proposed by 
Kawashima and Gubernatis pi] . Results from these two 
codes were cross-checked and agree for all spins within 
statistical errors, which gives a high degree of confidence 
to our calculations. 

The finite-volume correlation lengths S,{L) were cal- 
culated using a second-moment method, similar to 
Eq. (4.13) in Ref. |2|. Afterwards, the CEFPS finite- 
size-scaling method was applied. This method expresses 
S,{2L) / ^{L) as a universal function F{^{L)/ L) - appli- 
cable to all models within the universality class of the 
2-d lattice-regularized nearest-neighbour 0(3) non-linear 
cr-model. Iteration of F(^(L)/L) yields ^ = ^(L -^ 00). 
Unlike the case for S = 1/2, we found that for S > 1/2 
and our level of precision (« 0.2% for £,{L)), it is not nec- 
essary to incorporate any correction for scaling violations, 
even for lattice sizes as small as L/a ~ 10. The maximum 
^/a generated in our study are approximately 170,000 for 
S" = 1 and 135,000 for S = 5/2 [1^. Note that it is the 
finite-size scaling technique that enables the estimation 
of correlation lengths much larger than direct measure- 
ment allows (cf. Ref. pi] ], which achieved a maximum 
^/a = 24.94(7) at t = 0.18 on a square lattice with side 
length L/a — 200). Details of our algorithms, improved 
estimators, finite-size-scaling technique, and calculation 
of ^{L) will be provided in a separate paper p6|. 



Figure |l] shows the QMC data plotted on a Memphis 
chart, where we divide ^ by the leading (2-loop) term 
(e/8)(c/27rps)exp(l/i) and plot versus t. For S — 1/2, 
the agreement between QMC and CH3N2B down to 
ila « 10 (i Ri 0.3) is striking. For S > 1/2 we find that 
the QMC data smoothly merge into the CH3N2B predic- 
tions at C/a > 100 (t < 0.15) for S* = 1, and ^/a > 500 
{t < 0.10) for S = 5/2. The agreement in each case de- 
grades above some temperature. This is to be expected 
because CH3N2B leaves out higher-order terms from the 
CPT and spin- wave expansions. In addition, at high tem- 
peratures the field-theoretical requirements ^ :s> c/T and 
^ ^ a are no longer satisfied, and predictions such as 
CH3N2B become meaningless. In our figures we plot the 
CH3N2B predictions only for ^/a > 3. 

Ref. ^] found that for S = 1/2, the true value of spin 
stiffness ps is about 3% higher than the value predicted 
by third-order spin- wave theory (SW3). We find that 
comparison of QMC data with CH3N2B for S = 1/2 is 
in agreement with the fitted values of ps = 0.1800(5) 
and c = 1.657(2) found in Ref. |2^ (we set J and a to 
unity). That study combined the correlation length data 
fit to CH2N2, E g. (P), with a fit of finite- volume magnetic 
susceptibilities [ |l8yi9| to the predictions of CPT for the 
finite-size and temperature effects in the AFHM [^ . High 
fit precision was achieved by exploiting this combination. 
We will present a similar study for the full range of spins 
S < 5/2 in a separate paper p^ . 

For this Letter, we have choosen a different approach 
which only involves the correlation length. We demon- 
strate that for S > 1 one can directly rely on the SW3 
results [Oj to achieve a consistent connection between 
the QMC data and CH3N2B. For the 5 = 1 case 
we find that the SW3 predictions pf"'' = 0.869 and 
j^sw3 _ 3_Qg7 aj-g nearly correct. Our two-parameter fit 
gives pjpf^^ = 1.005(3) and c/c^^-' = 0.98(2). These ra- 

The 



tios correspond to ps = 0.8733(23) and c = 3.01(6) 
fit includes the ^/a > 100 {t < 0.15) data in Figure yj, and 
has x^/d.o.f. — 1.085 with 58 degrees of freedom (which 
corresponds to a significance level p = 30.5%). These 
values of ps and c are used in the figures for 5=1. 

Although the relative deviations from SW3 (0.5(3)% 
for Ps and —2(2)% for c) seem small, there is in fact a 
serious discrepancy. Using the SW3 values {pf"'-^ , c^™'-^) 
in CII3N2B and comparing to the QMC ^ gives x^/60 — 
1.39 (a poor fit, with p = 2.5%). Compared to the two- 
parameter fit, SW3 has Ax^ = -f20.2 (i.e., outside the 
99.99% confidence region). The reason the near-overlap 
with SW3 is deceptive is that the fit parameters ps and 
c are strongly anticorrelated (with correlation coefficient 
r — —0.977). Interestingly, the major axis of the nearly 
degenerate error ellipse is almost orthogonal to the curves 
of constant 9 = Ps/c^- In particular, the 68.3% confi- 
dence region for the joint probability distribution of ps 
and c (enclosed by the Ax^ = -1-2.30 ellipse) intersects 
the curve Ps/c^ = 0^""^. In other words, this 68.3% confi- 



1=1, we took from 



[|l) the 4.4% deviation between O""^'^/^ = 0.06277 



dence region contains {ps , c) values which are consistent 
with 6 = 6^^^ . We are thus motivated to check how close 
to pf"^ and c^™'' the corresponding single-parameter-fit 
values could in fact be. 

To do this, we set 6*5=1 = O^Zl = 0.09238 and per- 
formed a fit with the same set of 5 = 1 QMC data, with 
{ps, c) constrained to the one-dimensional parameter sub- 
space Ps = O^^^c^. As an upper bound for the error as- 
sociated with the assumption ^5=1 — 
Ref, ^ 

and that study's result ^5=1/2 = 0.06556. This choice 
is conservative since SWE is an expansion in powers of 
1/5, and is expected to become more accurate as spin in- 
creases. Upon refitting, we found ps/pf"^ = 1.0024(27), 
c/^swa ^ 1.001(21), and xV59 = 1.083 {p = 30.8%). 
These ratios correspond to ps = 0.8711(24) and c = 
3.07(6). (The difference between the one- and two- 
parameter fits would not be visible in Figure ^) Note 
the errors here are dominated by the conservative 4.4% 
uncertainty in 9s=i', the actual errors are bound to be 
smaller. 

For S = 5/2, we could not identify any deviation from 
the SW3 values pf"'' = 5.9444 and c^^^ = 7.3005. We 
found xVl4 = 1-194 {p = 27.2%) for the data with ^/a > 
500. These SW3 values of ps and c are used in the figures 
for S = 5/2. 

Figure shows the situation for S = 5/2 in more de- 
tail. There is an intermediate regime between ^/a w 500 
{t w 0.10), where CH3N2B starts to deviate, and £,/a w 
12 (t « 0.15), where the high-temperature-series expan- 
sion (HTE) llfl] starts to fail. Most of the experimental 
data on Rb2MnF4 fall into this "gap", which exists 
similarly for 5 = 1. At least for large spin S — 5/2, this 
intermediate regime is correctly described by the semi- 
classical approximation known as the pure-quantum self- 
consistent harmonic approximation (PQSCHA) |lj]. The 
diverse approaches collectively describe the S = 5/2 cor- 
relation length from extremely small to extremely large 
values. 

We note that a residual discrepancy between CII3N2B 
and numerical data persists in the classical limit S* — > 00, 
where the AFHM becomes the 2-d lattice-regularized 
nearest-neighbour 0(3) non-linear cr-model. Predictions 
for this model are available from analytical calculations 
H (^/a > 10^), Monte Carlo simulation |2|,||l (10 < 
^/a < 10^), and series expansion [0(^/a < 10). Hasen- 
fratz |l^] also supplies the 7—^0 form for the cor- 
rection exp(— (7(7)), which enables the computation of 
the classical S ^ 00 limit of Eq. (0). In Figure a, 
we plot the ratio of ^ to the CH3N2B prediction for 
S = 1/2,1,5/2, and cx) versus l/logio(0- By 5 = 5/2 
the discrepancy between the numerical data and CH3N2B 
has essentially reached the classical S ^ 00 limit. This 
means that the reasons for the residual discrepancy, in- 
cluding finite-order effects of the CPT and spin-wave ex- 
pansions, are the same for the quantum AFHM and the 



classical cr-model. 

In conclusion, the cutoff correction accounts for the 
overall spin dependence of the correlation length. The 
spin stiffness and spin-wave velocity approach the spin- 
wave theory predictions rapidly. The diverse approaches 
to the AFHM for higher-spin are complementary. 
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FIG. 1. Memphis chart for the AFHM. Correlation length 
QMC data and theoretical predictions are normalized by 
the leading (2-loop) prediction (e/8)(c/27rpi,) exp(l/f). Solid 
lines are the CH3N2B predictions, Eq. (kI). Dashed line is 
the spin- independent "no cutoff effects" prediction CH2N2, 
Eq. @. QMC data for S = 1/2 (diamonds), 5 = 1 (trian- 
gles) , and S — 5/2 (circles) ; solid (open) symbols are from the 
continuous (discrete) Euclidean-time algorithm. The QMC 
data smoothly merge into CH3N2B at low temperatures. 
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FIG. 2. Memphis chart for S = 5/2, showing CH3N2B 
(solid line); QMC (this work, filled circles); PQSCHA (dashed 
Une); high-temperature expansion (dotted line); and experi- 
ment (open squares). The QMC data connect the regimes 
where the various methods apply. 



